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We investigate the ground state properties of a 
two-dimensional electron gas in the lowest Lan- 
dau level using the Density Matrix Renormalization 
Group. The electron gas is confined to a cylinder 
with a strong magnetic field perpendicular to the 
surface. For a thin cylinder, the ground state is for 
generic filling factors a charge density wave state, 
however, atv = 1/2 a homogeneous and apparently 
gapless ground state is found and particle and hole 
excitations are studied. 

Of the many different methods used for study- 
ing strongly correlated systems the Density Ma- 
trix Renormalization Group (DMRG) invented by 
White |2j in 1992 is one of the most success- 
ful. The method, which is an extension of Wilson's 
Renormalization Group (RG) pj, provides a reli- 
able method for analyzing one-dimensional quan- 
tum systems that are far out of reach with exact 
diagonalization methods. For higher dimensional 
systems the DMRG is much less useful. 

The two-dimensional electron gas in a perpen- 
dicular magnetic field becomes effectively one- 
dimensional when the field is strong enough for the 
motion to be restricted to the lowest Landau level 
(LLL). In this limit the DMRG may become an 
efficient tool. The finite size DMRG has recently 
been applied to this system for an electron gas con- 
fined to a torus j^I- The results are in excellent 
agreement with exact diagonalization studies, ex- 
tending the latter from around 15 to 25 particles. 
In this article, we apply the infinite size version 
of the DMRG to the lowest Landau level problem 
when the electron gas is confined to a cylinder. This 
makes it possible to study the limit where the num- 



ber of particles N tends to infinity. On the other 
hand, the range of the interaction limits us to study 
thin cylinders. 

For generic quantum Hall (QH) filling fractions 
such asv = 1/3 we find charge density wave ground 
states where, for f — 1/3, approximately every 
third site is occupied by an electron. This result 
is in agreement with previous results on smaller 
systems [S] where the 1/3 state is shown to be a 
charge density wave on a thin cylinder. However, 
for the half-filled Landau level we find a non-trivial 
result — a homogeneous ground state with excita- 
tions that appear to be gapless. (For even smaller 
L, this state also gives way to a crystalline state.) 

We consider a charged particle moving in two 
dimensions in a perpendicular magnetic field B. 
By solving the Schrodinger equation in Landau 
gauge, A — Bxy, we map the two-dimensional sys- 
tem onto a one-dimensional system. The eigen- 
functions in the lowest Landau level are V'fe(^) — 



where L is the circumfer- 



ence of the cylinder, I — y |§ is the magnetic 
length and k is the momentum in the y-direction. 
Imposing periodic boundary conditions in the y- 
direction gives k = 2^ , m e Z and the den- 
sity of states in the LLL is = l/27r^^. The 
eigenfunction with j/-momentum k is centered at 
X = Xk = —kP. This provides a mapping of the 
LLL onto a one-dimensional problem which can be 
thought of as a spin chain where each site k has two 
states — either the site is occupied by an electron or 
it is empty. 

Following Trugman and Kivelson 0, we intro- 
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duce a repulsive short range potential between the 
particles which can be expanded as 

y(f)=^c,62^V2^52(f), (1) 

3 

where b is the range of the interaction. Taking the 
limit of short range, 6 — > 0, the expression (QJ re- 
duces to V{f) = ci6^V^(5^(r) since {S'^{f) ) = for 
any antisymmetric state. The form of this interac- 
tion can be motivated by the fact that the physics 
of quantum Hall systems does not change qualita- 
tively if the interaction is made more local than the 
more realistic Coulomb interaction. 

From the interaction and the single particle 
states in the LLL we find a second quantized form 
of the Hamiltonian. The Hamiltonian is then ana- 
lyzed and the low energy states are approximated 
by use of the DMRG method. The Hamiltonian is 

n k>l 

(2) 

where n, k and I denote all possible momentum 
states and we use units where £ — 1 and set = 1 . 
The DMRG algorithm becomes inefficient when the 
interaction extends over too many sites. The most 
local interaction in position space leads to the mo- 
mentum space interaction in We introduce a 
cut-off d in the range of the interaction — d is the 
number of sites the interaction extends over. Since 
the distance between the sites is 27r^^/L it follows 
that we are restricted to small circumferences L in 
order for the cut-off not to modify the shape of the 
interaction. Thus we are studying the QH system 
on a thin cylinder. Note, however, that alterna- 
tively we may view this as a large two-dimensional 
system with an interaction that extends only a dis- 
tance ~ 2TTdP/L in the a;-direction. Albeit being 
an unphysical interaction this is a mathematically 
well-defined limit where the system can be stud- 
ied. We perform calculations for d < 6. From the 
form of the interaction in (2) it follows that this 
restricts L to at most about 10 magnetic lengths. 
The results presented below are for L = 2n£ — this 
corresponds to the distance between the sites being 

e. 

As it stands, the Hamiltonian has a trivial 
ground state since the interaction is repulsive - the 
energy is minimal (and equal to zero) when there 



are no particles present. However, in a real physi- 
cal system the density of electrons. He, is finite — the 
negative charge of the electron gas is compensated 
for by a positive background charge that makes the 
whole system neutral. Thus we want to analyze 
the case when the filling fraction v = — = 2T:Pne 
is fixed and nonzero. In order to target a state 
with a specific filling fraction, v, we add a term to 
the Hamiltonian that favours the chosen density: 
H' = H + AH = H + ^N\j^YJ^^, cU - yf. (All 
energies presented below are of course for 7 = 0.) 

The Hamiltonian commutes with cjci, hence 
the total number of particles is a conserved quan- 
tum number. As a consequence, for given N 
only the filling fractions = 0, 1/iV, 2/7V, . . . , (iV - 
1)/A^, 1 can be obtained. This makes it in gen- 
eral necessary to keep states with different number 
of particles when studying a certain state in the 
N ^ 00 limit. 

For generic filling fractions we find charge den- 
sity wave ground states in agreement with earlier 
work As an example, we show in Fig. re- 

sults for ly — 1/3. The ground state is, to a good 
approximation, obtained by putting an electron on 
every third site. Thus the range of the interaction 
is so short that this simple crystal has lower energy 
than, e.g., a homogeneous state pj. 

However, for the half-filled Landau level we find a 
very different and much more interesting result — a 
homogenous ground state In Fig. 12 we present 
some of our results for the ;^ = 1/2 state. The ob- 
tained ground state has an almost perfectly uniform 
density throughout the system^ . The mean density 
is 1/2 up to the precision of at least 16 digits — this 
precision is a consequence of being a conserved 
quantum number. Furthermore, the standard de- 
viation of the single site density (cjci ) is as small 
as 10~^^ for large systems. We also note that the 
energy of the homogeneous state is very low com- 
pared to the crystalline states obtained when d <2. 
Thus the obtained state has used the off-diagonal 
elements of the Hamiltonian to reduce the energy 
in a very efficient way. The results quoted are inde- 
pendent of the number of states, m, that are kept 
in each iteration. 

The density ( cja ) given above, is really the lin- 

^This is true except near the boundaries of the system. 
The repulsive interaction, in combination with free boundary 
conditions, favours crystalline states for small systems and 
also near the edges of large systems. 
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Figure 1: The average energy per site (above) and the 
density distribution (cjci ) for 34 sites around the cen- 
ter of a large system (below) for the v = 1/3 ground 
state. 
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Figure 2: The f = 1/2 ground state energy per site 
(above) and the density distribution {c\ci) around the 
center of a large system (below). The density distri- 
bution inside the system is totally uniform. Further- 
more, the energy converges smoothly. That the energy 
increases is due to the fact that the edge effects become 
less important when the system size increases. 

ear density along the cylinder, i.e. the full density 
integrated around the cylinder. The density oper- 
ator is 

m,n 



m,n 

However, we find that (c]„c„) oc 6m, n to a very 
high precision^. Hence the density is constant on 
the cylinder^. 

For d < 2, the i/ — 1/2 ground state is crystalline. 
At c? = 1 every other site is occupied and at c? = 2 
two occupied sites alternate with two empty — these 
are exact ground states. At f — 1/4 the ground 
state is a crystal for d <6. 

We have made a preliminary investigation of 
charged excitations at = 1/2 by targeting states 
with density ly = 1/2±1/N for N sites^. We find a 
particle excitation with charge — |e| relative to the 
ground state and energy Ep — 2.0 and a hole with 
charge |e| and energy Eh = —1.8 {J = 1, d = A 
and we keep to = 4 states in total). This gives 
an energy gap A = Ep + Eh — 0.2 to creating a 
separated particle-hole pair. We take the smallness 
of A compared to Ep , Eh to indicate that A ap- 
proaches zero in the thermodynamic limit and that 
the system is gapless. To test this we should con- 
sider the limit N ^ oo and m — > oo. For the par- 
ticle this is straightforward, however, for the hole 
the calculations are unstable when m is increased. 
Increasing N only moderately (to avoid the insta- 
bility) we find that A decreases with increasing m, 
supporting the conclusion that A ^ 0. We believe 
the instability is due to a degeneracy in the energy 
of the hole state. We find that the hole forms at the 
center of the system — there should then be degen- 
erate states where the hole is translated to other 
positions. That no instability is seen for the parti- 
cle indicates that it forms at the edge of the system. 
Fig. El shows the density profile of the hole. The 
charge integrated over the area, with length 34£, 
shown in the figure is 0.98|e| relative to the ground 
state in Fig. [2 

A further indication of the existence of gapless 
excitations is shown in Fig. ^ It shows the in- 
crease in total energy of the v = 1/2 ground state 
per iteration up to large systems'^. In these calcula- 
tions, states at density i/ =1/2 and v = l/2-|-l/A^ 

^The deviation from this result is roughly of the same 
order as the deviation of the single site density (cjci), i.e. 
10-12. 

^The finite distance 2-kI/? jL = I between two sites leads 
to a variation in the spatial density of roughly 0.01%. 
■^When doing this we target also the = 1/2 state. 
'The figure shows only data from every second iteration 
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Figure 3: The density profile of a charge one hole in 
the V = 1/2 state. The plot shows the expectation values 
of the single site density operators (cjci ) for a range of 
sites near the center of the superblock. The line is a 
guide to the eye. 



Figure 4: Energy added to ground state in each itera- 
tion for half -filled LLL keeping states with v = 1/2 and 
V — 1/2 + 1/N . Initially, the ground state is homo- 
geneous, then eventually a transition to a state with a 
well separated particle-hole pair is observed. 



are kept. Initially, a smooth convergence to the ho- 
mogeneous V = \/2 ground state is observed — this 
is reflected in the flgure in that the same amount 
of energy is added at each step. In each iteration 
two sites and one electron are added at the center 
of the superblock, thus keeping the average density 
flxed at V = 1/2. Then suddenly at iteration 470 
a slightly smaller energy is added. A new ground 
state is formed as a superposition of the homoge- 
neous u = \/2 ground state in the previous itera- 
tion (adding two sites and one electron to this) and 
the state containing one particle excitation (adding 
two sites and no particle to this). The energy of 
this state is actually lower than the energy of the 
homogeneous state since the added energy in the 
last iteration decreased^. We take this as indicat- 
ing that the particle-hole excitation is gapless. 

In addition to the above studies of excitations, 
where we always target only the lowest energy 
state, we have also targeted the two or three low- 
est states. We then flnd that these states become 
degenerate as N and to ^ oo thus lending further 
support to the gaplessness of the ground state. 

since there is a "parity-effect" that distinguishes even and 
odd iterations. 

^Eventually, after further iterations, the state converges 
to a state with a well separated particle-hole pair. That this 
is picked out rather than the homogeneous state (or a linear 
combination) is a consequence of our algorithm. 



To conclude, for a thin cy Under we find a ho- 
mogeneous, presumably gapless, ground state for 
the half-filled lowest Landau level — at all other fill- 
ing fractions that we have investigated, including 

= 1/4, we find charge density wave states. It is 
of course tempting to identify the obtained ground 
state as the Rezayi-Read state that is known 
from comparison to exact diagonalization studies to 
describe the v = 1/2 state on the sphere very well. 
It remains to compare our state to the Rezayi-Read 
state. It is computationally nontrivial to compare 
large systems whereas small systems suffer from 
large edge effects, presumably making a compar- 
ison meaningless. A natural solution to this would 
be to employ the finite size DMRG method, possi- 
bly on the torus. An interesting questions is why 
there is a homogeneous state on thin cylinders only, 
as it seems, at v = 1/2. If this state is indeed the 
Rezayi-Read state, does our results shed any fur- 
ther light on the i/ ~ 1/2 quantum Hall state? A 
study of our i/ = 1/2 ground state indicates that it 
can be understood in terms of free fermions [TT] . 
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